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Abstract. Type la supernovae (SN la) are generally believed to be the result of the 
thermonuclear disruption of Chandrasekhar-mass carbon-oxygen white dwarfs, mainly 
because such thermonuclear explosions can account for the right amount of 56 Ni, which 
is needed to explain the light curves and the late-time spectra, and the abundances of 
intermediate-mass nuclei which dominate the spectra near maximum light. Because of 
■^j- ' their enormous brightness and apparent homogeneity SN la have become an important 

tool to measure cosmological parameters. 

In this article the present understanding of the physics of thermonuclear explosions 
is reviewed. In particular, we focus our attention on subsonic ("deflagration") fronts, 
| i.e. we investigate fronts propagating by heat diffusion and convection rather than by 

compression. Models based upon this mode of nuclear burning have been applied very 
successfully to the SN la problem, and are able to reproduce many of their observed 
features remarkably well. However, the models also indicate that SN la may differ 
considerably from each other, which is of importance if they are to be used as standard 
candles. 
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INTRODUCTION 

Type la supernovae, i.e. stellar explosions which do not have hydrogen in their 
spectra, but intermediate-mass elements, such as silicon, calcium, cobalt, and iron, 
have recently received considerable attention because it appears that they can be 
used as "standard candles" to measure cosmic distances out to billions of light 
years away from us. Moreover, observations of type la supernovae seem to indicate 
that we are living in a universe that started to accelerate its expansion when it was 
about half its present age. These conclusions rest primarily on phenomenological 
models which, however, lack proper theoretical understanding, mainly because the 
explosion process, initiated by thermonuclear fusion of carbon and oxygen into 
heavier elements, are difficult to simulate even on supercomputers, for reasons we 
shall discuss in this article. 

The most popular progenitor model for the average type la supernovae is a 
massive white dwarf, consisting of carbon and oxygen, which approaches the Chan- 



drasekhar mass, M Ch ~ 1.39 M Q , by a yet unknown mechanism, presumably 
accretion from a companion star, and is disrupted by a thermonuclear explosion 
(see, e.g., [56] for a review). Arguments in favor of this hypothesis include the 
ability of the models to fit the observed spectra and light curves rather well. 

However, not only is the evolution of massive white dwarfs to explosion very 
uncertain, leaving room for some diversity in the allowed set of initial conditions 
(such as the temperature profile at ignition), but also the physics of thermonuclear 
burning in degenerate matter is complex and not well understood. The generally 
accepted scenario is that explosive carbon burning is ignited either at the center of 
the star or off-center in a couple of ignition spots, depending on the details of the 
previous evolution. After ignition, the flame is thought to propagate through the 
star as a sub-sonic deflagration wave which may or may not change into a detonation 
at low densities (around 10 7 g/cm 3 ). Numerical models with parameterized velocity 
of the burning front have been very successful, the prototype being the W7 model 
of [37]. However, these models do not solve the problem because attempts to 
determine the effective flame velocity from direct numerical simulations failed and 
gave velocities far too low for successful explosions [25,18,1]. This has led to some 
speculations about ways to change the deflagration into a supersonic detonation 
[14,15]. All these aspects of supernova models will be discussed in the following 
sections. 

NUCLEAR BURNING IN DEGENERATE C+O 

MATTER 

Owing to the strong temperature dependence of the nuclear reaction rates nu- 
clear burning during the explosion is confined to microscopically thin layers that 
propagate either conductively as subsonic deflagrations ( "flames" ) or by shock com- 
pression as supersonic detonations [6,21]. Both modes are hydrodynamically un- 
stable to spatial perturbations as can be shown by linear perturbation analysis. In 
the nonlinear regime, the burning fronts are either stabilized by forming a cellular 
structure or become fully turbulent - either way, the total burning rate increases 
as a result of flame surface growth [24,53,60]. Neither flames nor detonations can 
be resolved in explosion simulations on stellar scales and therefore have to be rep- 
resented by numerical models. 

When the fuel exceeds a critical temperature T c where burning proceeds nearly 
instantaneously compared with the fluid motions a thin reaction zone forms at 
the interface between burned and unburned material. It propagates into the sur- 
rounding fuel by one of two mechanisms allowed by the Rankine-Hugoniot jump 
conditions: a deflagration ("flame") or a detonation. 

If the overpressure created by the heat of the burning products is sufficiently 
high, a hydrodynamical shock wave forms that ignites the fuel by compressional 
heating. Such a self-sustaining combustion front that propagates by shock-heating 
is called a detonation. If, on the other hand, the initial overpressure is too weak, 



the temperature gradient at the fuel-ashes interface steepens until an equilibrium 
between heat diffusion and energy generation is reached. The resulting combustion 
front consists of a diffusion zone that heats up the fuel to T c , followed by a thin 
reaction layer where the fuel is consumed and energy is generated. It is called a 
deflagration or simply a flame and moves subsonically with respect to the unburned 
material [21]. Flames, unlike detonations, may therefore be strongly affected by 
turbulent velocity fluctuations of the fuel. Only if the unburned material is at 
rest, a unique laminar flame speed S\ can be found which depends on the detailed 
interaction of burning and diffusion within the flame region, e.g. [60]. Following 
[21], it can be estimated by assuming that in order for burning and diffusion to 
be in equilibrium, the respective time scales, Tb ~ e/w and ~ S 2 /k, where S 
is the flame thickness and k is the thermal diffusivity, must be similar: Tb ~ 
Defining S\ = one finds S\ ~ (nw/e) 1 ^ 2 , where w should be evaluated at 

T T c [50]. This is only a crude estimate due to the strong T-dependence of 
w. Numerical solutions of the full equations of hydrodynamics including nuclear 
energy generation and heat diffusion are needed to obtain more accurate values 
for S\ as a function of p and fuel composition. Laminar thermonuclear carbon 
and oxygen flames at high to intermediate densities were investigated by [4,12,54], 
and, using a variety of different techniques and nuclear networks, by [50]. For the 
purpose of SN la explosion modeling, one needs to know the laminar flame speed 
Si ~ 10 7 . . . 10 4 cm s _1 for p ss 10 9 . . . 10 7 g cm -3 , the flame thickness 5 = 10~ 4 ... 1 
cm (defined here as the width of the thermal pre-heating layer ahead of the much 
thinner reaction front), and the density contrast between burned and unburned 
material p = Ap/p = 0.2. . .0.5 (all values quoted here assume a composition of 
X c = X Q = 0.5 [50]). The thermal expansion parameter p reflects the partial 
lifting of electron degeneracy in the burning products, and is much lower than the 
typical value found in chemical, ideal gas systems [53]. 

HYDRODYNAMIC INSTABILITIES AND 
TURBULENCE 

The best studied and probably most important hydrodynamical effect for mod- 
eling SN la explosions is the Rayleigh- Taylor (RT) instability resulting from the 
buoyancy of hot, burned fluid with respect to the dense, unburned material 
[29,30,25,17,18,33], and after more than five decades of experimental and numerical 
work, the basic phenomenology of nonlinear RT mixing is fairly well understood 
[8,22,47,43,57]: Subject to the RT instability, small surface perturbations grow un- 
til they form bubbles (or "mushrooms") that begin to float upward while spikes 
of dense fluid fall down. In the nonlinear regime, bubbles of various sizes interact 
and create a foamy RT mixing layer whose vertical extent /irt grows with time t 
according to a self-similar growth law, h RT = ag(p/2)t 2 , where a is a dimensionless 
constant (a ~ 0.05) and g is the background gravitational acceleration. 

Secondary instabilities related to the velocity shear along the bubble surfaces [34] 



quickly lead to the production of turbulent velocity fluctuations that cascade from 
the size of the largest bubbles (« 10 7 cm) down to the microscopic Kolmogorov 
scale, Zk R3 10~ 4 cm where they are dissipated [33,18]. Since no computer is capable 
of resolving this range of scales, one has to resort to statistical or scaling approxi- 
mations of those length scales that are not properly resolved. The most prominent 
scaling relation in turbulence research is Kolmogorov's law for the cascade of ve- 
locity fluctuations, stating that in the case of isotropy and statistical stationarity, 
the mean velocity v of turbulent eddies with size / scales as v ~ Z 1//3 [20]. Given 
the velocity of large eddies, e.g. from computer simulations, one can use this rela- 
tion to extrapolate the eddy velocity distribution down to smaller scales under the 
assumption of isotropic, fully developed turbulence [33]. Knowledge of the eddy 
velocity as a function of length scale is important to classify the burning regime 
of the turbulent combustion front [36,35,19]. The ratio of the laminar flame speed 
and the turbulent velocity on the scale of the flame thickness, K = S\/v(S), plays 
an important role: if K >> 1, the laminar flame structure is nearly unaffected by 
turbulent fluctuations. Turbulence does, however, wrinkle and deform the flame on 
scales / where Si <C v(l), i.e. above the Gibson scale l g defined by S\ = v(l g ) [40]. 
These wrinkles increase the flame surface area and therefore the total energy gen- 
eration rate of the turbulent front [7] . In other words, the turbulent flame speed, 
S t , defined as the mean overall propagation velocity of the turbulent flame front, 
becomes larger than the laminar speed S\. If the turbulence is sufficiently strong, 
v(L) 3> Si, the turbulent flame speed becomes independent of the laminar speed, 
and therefore of the microphysics of burning and diffusion, and scales only with the 
velocity of the largest turbulent eddy [7,5]: 



Because of the unperturbed laminar flame properties on very small scales, and the 
wrinkling of the flame on large scales, the burning regime where K S> 1 is called 
the corrugated flamelet regime [42,5]. 

As the density of the white dwarf material declines and the laminar flamelets be- 
come slower and thicker, it is plausible that at some point turbulence significantly 
alters the thermal flame structure [19,36]. This marks the end of the flamelet 
regime and the beginning of the distributed burning, or distributed reaction zone, 
regime, e.g. [42]. So far, modeling the distributed burning regime in exploding 
white dwarfs has not been attempted explicitly since neither nuclear burning and 
diffusion nor turbulent mixing can be properly described by simplified prescriptions 
(see, however, [27]). Phenomenologically, the laminar flame structure is believed 
to be disrupted by turbulence and to form a distribution of reaction zones with 
various lengths and thicknesses. In order to find the critical density for the tran- 
sition between both regimes, we need to formulate a specific criterion for flamelet 
breakdown. A criterion for the transition between both regimes is discussed in 



S t ~ v(L) . 



(1) 



[36,35] and [19]: 



^cutoff < 8 ■ 



(2) 



Inserting the results of [50] for S\ and 5 as functions of density, and using a typ- 
ical turbulence velocity t>(10 6 cm) ~ 10 7 cm s -1 , the transition from flamelet to 
distributed burning can be shown to occur at a density of p dis « 10 7 g cm -3 [35]. 

MODELING TURBULENT THERMONUCLEAR 

COMBUSTION 

Next we will outline a way by which several of the ideas discussed in the previ- 
ous sections can be incorporated into a numerical scheme to model thermonuclear 
combustion in SN la. 

Numerical simulations of any kind of turbulent combustion have always been a 
challenge, mainly because of the large range of length scales involved. In type la 
supernovae, in particular, the length scales of relevant physical processes range from 
10 _3 cm for the Kolmogorov-scale to several 10 7 cm for typical convective motions. In 
the currently favored scenario the explosion starts as a deflagration near the center 
of the star. Rayleigh- Taylor unstable blobs of hot burnt material are thought to 
rise and to lead to shear-induced turbulence at their interface with the unburnt gas. 
This turbulence increases the effective surface area of the flamelets and, thereby, 
the rate of fuel consumption; the hope is that finally a fast deflagration might 
result, in agreement with phenomeno logical models of type la explosions [37]. 

Despite considerable progress in the field of modeling turbulent combustion for 
astrophysical flows (see, e.g., [33]), the correct numerical representation of the ther- 
monuclear deflagration front is still a weakness of the simulations. Methods used 
up to now are based on for the reactive-diffusive flame model [16], which artificially 
stretches the burning region over several grid zones to ensure an isotropic flame 
propagation speed. However, the soft transition from fuel to ashes stabilizes the 
front against hydrodynamical instabilities on small length scales, which in turn re- 
sults in an underestimation of the flame surface area and - consequently - of the 
total energy generation rate. Moreover, because nuclear fusion rates depend on 
temperature nearly exponentially, one cannot use the zone-averaged values of the 
temperature obtained this way to calculate the reaction kinetics. 

Therefore we have decided to use a front tracking method to cure some of these 
weaknesses. The method is based on the so-called level set technique which was 
originally introduced by Osher and Sethian [38]. They used the zero level set of 
a n-dimensional scalar function to represent (n — l)-dimensional front geometries. 
Equations for the time evolution of such a level set which is passively advected by 
a flow field are given in [49]. The method has been extended to allow the track- 
ing of fronts propagating normal to themselves, e.g. deflagrations and detonations 
[48,45,44]. In contrast to the artificial broadening of the flame in the reaction- 
diffusion-approach, this algorithm is able to treat the front as an exact hydrody- 
namical discontinuity. We will demonstrate that such a method can be applied to 
the supernova problem and, in addition, we will show that even if one attempts 
to model the physics of thermonuclear burning on unresolved length scales well 
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5 Pf©URE 1. Snapshots of the temperature and the front geometry at 0.4s for model B5 of 
Reinecke et al. (1999a) (right figure) and turbulent velocity fluctuations on the grid scale (left 
panel). The position of of the front is indicated by the dotted curve. 

by physically motivated "Large Eddy Simulations" (LES), one still has to perform 
calculations with very high spatial resolution in order to be at least near to a con- 
verged solution. The numerical results presented here were obtained by performing 
2D simulations, but there is no problem extending them to 3D. Such simulations 
are under way. 

APPLICATION TO THE SUPERNOVA PROBLEM 

We have carried out numerical simulations in cylindrical rather than in spherical 
coordinates, mainly because it is much simpler to implement the level set on a 
Cartesian (r,z) grid. Moreover, the CFL condition is somewhat relaxed in compari- 
son to spherical coordinates. The grid we used in most of our simulations maps the 
white dwarf onto 256x256 mesh points, equally spaced for the innermost 226x226 
zones by A=1.5-10 6 cm, but increasing by 10% from zone to zone in the outer parts. 
The white dwarf, constructed in hydrostatic equilibrium for a realistic equation of 
state, has a central density of 2.9-10 9 g/cm 3 , a radius of 1.5-10 8 cm, and a mass of 
2.8-10 33 g, identical to the one used by [33]. The initial mass fractions of C and O 
are chosen to be equal, and the total binding energy turns out to be 5.4-10 50 erg. 
At low densities (p < 10 7 g/cm 3 ), the burning velocity of the front is set equal 
to zero because the flame enters the distributed regime and our physical model is 
no longer valid. However, since in reality some matter may still burn the energy 
release obtained in the simulations is probably somewhat too low. 

Here we present only the results for two models, one in which nuclear burning 
was ignited off-center in a blob and a second one in which initially five blobs were 
burned as an initial condition, and refer to [44] for simulations with other initial 
conditions. 

First, in Fig.l a snapshot of the "five-blob" model B5 is shown at t = 0.4s. 
The left panel gives the position of the burning front, the temperatures (in gray- 
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FIGURE 2. Snapshots of the temperature and the front geometry at 1.05s, taken for the 
low-resolution model (left figure) and high-resolution run, respectively. 



shading), and the expansion velocities. The right panel shows the distribution of 
turbulent velocity fluctuations. We find that, in accord with one's intuition, most 
of the turbulence is generated in a very thin layer near the front. Since in the 
limit of high turbulence intensity the nuclear flames propagate with the turbulent 
velocity it is obvious that this propagation velocity exceeds the laminar flame speed. 
However, for most of the initial conditions we have investigated this increase was 
not sufficient to unbind the star. In fact, model B5 of [44] was the only one of the 
set that did explode. 

Moreover, we show the results obtained with 3 times higher resolution for com- 
parison. Fig. 2 gives a snapshot taken at 1.05s after ignition for the low- resolution 
run (left figure) and the high-resolution run, respectively. Although the increase 
in spatial resolution is only a factor of 3, the right panel shows clearly more struc- 
ture. This is an important effect because in the flamelet regime the rate of fuel 
consumption, in first order, increases proportional to surface area of the burning 
front. The net effect is that the low-resolution model stays bound at the end of 
the computations, whereas the better resolved model explodes with an explosion 
energy of about 2xl0 50 erg. Fig. 2 also demonstrates that the level-set prescrip- 
tion allows to resolve the structure of the burning front down almost to the grid 
scale, thus avoiding artificial smearing of the front which is an inherent problem 
of front-capturing schemes. We want to stress that this gain of accuracy is not 
obtained at the expense of smaller CFL time steps because in our hybrid scheme 
the hydrodynamics is still done with cell-averaged quantities. 

We also want to stress again, however, that the numerical simulations described 
here, although leaving behind unbound white dwarfs in several cases, would not 
fit the observed spectra of SN la, mainly because the expansion velocity of the 
partially burnt gas is too low. Whether this problem is cured once the model 
calculations are extended to three spatial dimensions has to be seen, but one can 
expect significant changes because of the different ratios of volume to surface area 
of burning blobs in 3D. 



DELAYED DETONATIONS? 



Turbulent deflagrations can sometimes be observed to undergo spontaneous tran- 
sitions to detonations (deflagration-detonation transitions, DDTs) in terrestrial 
combustion experiments [53]. Thus inspired, it was suggested that DDTs may 
occur in the late phase of the explosion, providing an elegant explanation for the 
initial slow burning required to pre-expand the star, followed by a fast combustion 
mode that produces large amounts of high-velocity intermediate mass elements 
[14,55]. Many ID simulations have meanwhile demonstrated the capability of the 
delayed detonation scenario to provide excellent fits to SN la spectra and light 
curves [56,11], as well as reasonable nucleosynthesis products with regard to solar 
abundances [15,13]. In the best fit models, the initial flame phase has a rather slow 
velocity of roughly one percent of the sound speed and transitions to detonation at 
a density of Pddt ~ 10 7 g cm' 3 [11,13]. The transition density was also found to 
be a convenient parameter to explain the observed sequence of explosion strengths 

[11]- 

Various mechanisms for DDT were discussed in the early literature on delayed 
detonations (see [36] and references therein). Recent investigations have focussed 
on the induction time gradient mechanism [59,23], analyzed in the context of SNe 
la by [2] and [3]. It was realized by [19] and [36] that a necessary criterion for this 
mechanism is the local disruption of the flame sheet by turbulent eddies, or, in 
other words, the transition of the burning regime from "flamelet" to "distributed" 
burning. Simple estimates [35] show that this transition should occur at roughly 
10 7 g cm -3 , providing a plausible explanation for the delay of the detonation. 

The critical length (or mass) scale over which the temperature gradient must 
be held fixed in order to allow the spontaneous combustion wave to turn into a 
detonation was computed by [19] and [36]; it is a few orders of magnitude thicker 
than the final detonation front and depends very sensitively on composition and 
density. Recently, this conclusion was confirmed by [27], who computed the in- 
teraction between turbulence and the nuclear flame in the distributed regime by 
means of a particular turbulence model, and by [26]. 

The virtues of the delayed detonation scenario can be summarized as follows. It is 
undisputed that suitably tuned delayed detonations satisfy all the constraints given 
by SN la spectra, light curves, and nucleosynthesis. If Pddt is indeed determined 
by the transition of burning regimes - which in turn might be composition depen- 
dent [52] - the scenario is also fairly robust and Pddt m &y represent the explosion 
strength parameter. Note that in this case, the variability induced by multi-point 
ignition needs to be explained away. If, on the other hand, thermonuclear flames 
are confirmed to be almost unquenchable, the favorite mechanism for DDTs be- 
comes questionable [31]. Moreover, should the mechanism DDT rely on rare, strong 
turbulent fluctuations one must ask about those events that fail to ignite a deto- 
nation following the slow deflagration phase which, on its own, cannot give rise to 
a viable SN la explosion. They might end up as pulsational delayed detonations or 
as unobservably dim, as yet unclassified explosions. Multidimensional simulations 



of the turbulent flame phase may soon answer whether the turbulent flame speed 
is closer to 1 % or 30 % of the speed of sound and hence decide whether DDTs are 
a necessary ingredient of SN la explosion models. 



SUMMARY AND CONCLUSIONS 



In this review we have outlined our present understanding of Type la supernovae 
explosions. From the tremendous amount of work carried out over the last couple 
of years it has become obvious that the physics of SNe la is very complex, ranging 
from the possibility of very different progenitors to the complexity of the physics 
leading to the explosion and the complicated processes which couple the interior 
physics to observable quantities. None of these problems is fully understood yet, 
but what one is tempted to state is that, from a theorist's point of view, it appears 
to be a miracle that all the complexity seems to average out in a mysterious way 
to make the class so homogeneous. In contrast, as it stands, a save prediction from 
theory seems to be that SNe la should get more divers with increasing observed 
sample sizes. If, however, homogeneity would continue to hold this would certainly 
add support to the Chandrasekhar-mass scenario discussed in this article. On the 
other hand, even an increasing diversity would not rule out Chandrasekhar-mass 
progenitors for most of them. In contrast, there are ways to explain how the 
diversity is absorbed in a one parameter family of transformations, such as the 
Phillips- relation [41,10] or modifications of it [46,39]. 

As far as the explosion/combustion physics and the numerical simulations are 
concerned significant recent progress has made the models more realistic (and re- 
liable). Thanks to ever increasing computer resources 3-dimensional simulations 
have become feasible which treat the full star with good spatial resolution and 
realistic input physics. Already the results of 2-dimensional simulations indicate 
that pure deflagrations waves in Chandrasekhar-mass C+O white dwarfs can lead 
to explosions, and one can expect that going to three dimensions, because of the 
increasing surface area of the nuclear flames, should add to the explosion energy. 
If confirmed, this would eliminate pulsational detonations from the list of potential 
models. On the side of the combustion physics, the burning in the distributed 
regime at low densities needs to be explored further, but it is not clear anymore 
whether a transition from a deflagration to a detonation in that regime is needed 
for successful models. In fact, according to recent studies such a transition appears 
to be rather unlikely. 
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